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Abstract 

We construct the phase diagrams of the quark-antiquark and diquark condensates at finite tem- 
perature and density in the 3D (dimensional) 2-flavor Gross Neveu model. We found that, in 
contrast to the case of the 4D Nambu Jona-Lasinio model, there is no region where the quark- 
antiquark and diquark condensates coexist. The phase diagrams obtained for some parameter 
region show similar structure with the 4D QCD phase diagram. 
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I. INTRODUCTION 



The four-fermion interaction models are considered to be the effective theory for describ- 
ing phase transitions. In 3 + 1 dimensions (4D), the Nambu Jona-Lasinio (NJL) model can 
be regarded as an effective theory of QCD l(. The model successfully reproduces the chiral 



phase transition in quark matter (see, e.g. |2|,|3|]). Furthermore, through recent studies of the 
diquark condensate, a variety of color superconducting phases are expected to be realized[4j. 
It has been shown that there is the region where the chiral (quark- antiquark) and diquark 
condensates coexist. (For a nice review, see, e.g. ^j.) 

The Gross-Neveu (GN) modelproposed in 1974, is a model of Dirac fermions interact- 
ing via four-fermion interactions jg]. The 2D GN model is a renormalizable quantum field 
theory, and the 3D GN model is renormalizable in the leading 1/N order. They are closely 
related to the Bardeen-Cooper-Schrieffer theory of superconductivity [7J. Although the lower 
dimensional theory seems not to be realistic, the 2D GN type models are believed to be the 
effective models of 1 dimensional condensed matter systems such as conducting polymers like 
polyacetylene[8]. Since the GN model shares many properties with QCD, notably asymp- 
totic freedom, chiral symmetry breaking in vacuum, the comparison of QCD and the GN 
model is an interesting subject. 



A variety of works has been devoted to the study of the GN model] 



0,Q,[iil,Q,Q,Q,H, 



Hi. By using bag-node, boundary conditions, a Cosed fornuda for the effective _li Z ed 
coupling constant in the large-N limit, was derived in [10]. In a constant curvature space, 
the phase structure of chiral symmetry breaking in the GN model at finite temperature 

n 

and density is discussed in [llj]. The phase diagrams of the quark- antiquark condensate in 



the 2D GN model was obtained in 



121 ] . In 3D, analyses so far have been done choosing 



2-dimensional(2d) or 4-dimensional(4d) spinor representations for quarks. The case of the 
2d representation is interesting itself because it is the nontrivial lowest-order representation. 
The phase diagram of the quark- antiquark condensate was constructed by employing the 2d 
representation i n [131] . On the other hand in the case of the 4d representation, there exists 
the 7 5 (see, e.g. |l7|) and the properties of the model bear resemblance to the NJL model in 
4D. In this sense, this case is also interesting and the phase diagram of the quark- antiquark 
condensate was obtained in [14J]. The relation between the 2d and 4d representations are 
discussed in [15]. Furthermore, through analyzing the 3D GN with the 2d representation, 
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the phase structure of the quark-antiquark(gg) and diquark(gg) condensates in vacuum (zero 
temperature and chemical potential(T = /x = 0)) was studied in [3] under the condition 
that the qq and qq condensates are much smaller than the cut-off scale of the model. It was 
found that the qq and qq condensates do not coexist at T = fi — 0. 

In this paper, we study the qq and qq condensates at finite temperature and density in 
the 3D GN model with the 2d representation. We obtain the phase diagrams, and discuss 
the similarities and differences between the 3D GN model and the 4D NJL model. 

The plan of the paper is as follows: In Sec. II we present the Lagrangian of the 3D GN 
model and introduce the mean-field approximation. In Sec. Ill we derive the thermodynamic 
potential. In Sec. IV we present the results of the numerical analyses. Among others, we 
show that there is no region where the qq and qq condensates coexist. In Sec.V we display the 
phase diagrams. We find that the structure of phase diagrams bear resemblance to the QCD 
phase diagram for some parameter region. Sec. VI is devoted to summary and conclusions. 
In Appendixes A and B we describe the intermediate calculation and the renormalization 
for the thermodynamic potential. 



II. GROSS NEVEU MODEL 

The general form of the Lagrangian density of the 3D 2 flavor massless Gross Neveu 
model with 2d representation reads 

C = qifiq + G s [(qq) 2 + (qrq) 2 ] 
8 

+ G D J2(Qr\aq C )(fr\ a q). (1) 

a=l 

Here f — (ti, T2, 73) are the Pauli matrices in flavor space and A a is a Gell-Mann matrix in 
color space. Gs and Gn are the coupling constants of the qq and qq interactions, respectively, 
and C is the charge conjugation matrix. For the Dirac 7 matrices and C, we use the forms 
as in . 

*-(::)"-(::)■*-(::)-* 

The charge conjugated field are defined by 

q c = C f, f = q T C. (3) 



Following the reasonings described in [5|, we reduce the Lagrangian density in Eq.flT]) to 

£ = qiftq + G s {qqf 

+ G D {qT 2 \ a q c ){fr 2 \ a q). (4) 

a=2,5,7 

C in Eq.(j4"]) enjoys various symmetry properties, which are fully discussed in [if]]. Here we 
choose a color direction for diquark condensate to blue, which is equivalent to select A2 in 
Eq.fllj) (see [ij]). Then we finally arrive at the following Lagrangian density: 

C = qiftq + G s {qq) 2 

+ G D (qT 2 X 2 q c )(q C T 2 X 2 q). (5) 

Due to A2 in Eq.Q, only two colors (red, green) participate in the qq condensate while all 
three colors (red, green, blue) do in the qq condensate. 

Let us introduce the mean-field approximation and rewrite the Lagrangian density as 
follows: 

C =qi$q - qaq - ^A*(q c r 2 X 2 q) 

1 A / , o 2 |A| 2 

where a and A are the order parameters for the qq and qq condensates: 

a = -2G s (qq) and A = -2G D (q c r 2 X 2 q). (7) 

To deal with finite density system, we introduce a chemical potential being conjugate to 
the quark number q^q. Then the Lagrangian density reads 

£ =qi$q + q/ij°q - qaq - ^A*(q c r 2 X 2 q) 



Introducing the Nambu-Gorkov basis 18| 



= j q I and * = ( q q c ) 

q c 

and using the relation q Cr y°q c = —q^q, we can write the Lagrangian density in a momentum 
space as 
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where 



G- 1 



(jf-<7 + /i7 )l/l c -r 2 A 2 Al s 
-r 2 A 2 A*l s {tf - a - w°)l f l c 



(10) 



If, l c and l s are the unit matrix in flavor, color and spinor space respectively. 



III. THE THERMODYNAMIC POTENTIAL 



A. Derivation of the thermodynamic potential 

Following the standard method, we can evaluate the thermodynamic potential: 

Q(a,\A\) ^ UM 



AGs AG 



D 



-Lin Ad*]exp i^C/KT 1 )* 

L n,p 



(11) 



where j3 = 1/T is the inverse temperature and V is the volume of the system. With the 
help of the formula 



[d^f] exp 

we can rewrite Eq. (TTT|) as 

Q(a, |A| 



Ad*] exp i^T^CT 1 )* 



a 



IAI 2 1 



Det^iPG- 1 ) 



lnDet 1 / 2 ^- 1 ). 



4G S 4G D PV 

After some manipulations which is given in Appendix A, the determinant becomes 

Det^iG- 1 ) = [p 2 -EX 2 ] 2 [pt-E- A 2 } 2 [pt-E^}[pt-E-% 

where po = i(2n + l)vrT, (n = • • • , —2, —1, 0, 1, 2, • • • ) and 

E ± = E±/i, E = 2 + a 2 , p 2 =p 2 +p 2 
£± 2 = £ 2 + /i 2 + |A| 2 ± 2 v / J EV + a 2 |A| 2 (> 0). 



Thus, we obtain 



n(a,\A\) 



IAI 



AG S AGd 
d 2 p 



7 S^./ (2-) 2 



ln[/3 2 (p 2 - J E; ±2 )]+21n[/? 2 (p 2 - J E;f)] 



(12) 



(13) 



(14) 



(15) 



(16) 



The frequency summation may be performed in a standard manner 191 ]: 



HP 2 (pI - E 2 )} = P[E + 2Tln(l + e^ E )\. 



(17) 



Then we finally obtain 



Q(a, | A|) = Q (a, |A|) + Q T (cr, |A 
- 2 |A| 2 r A2 



fio(t7,|A|) 



<7 



+ 



fi T ( C r,iAi)=-2T^ y. 



d 2 p 



hy(l + e- pE ) + 21n(l + e-^A; 



(18) 
(19) 

(20) 



Here is T independent contribution, which is ultraviolet divergent, while the temperature 
dependent part Qt is finite. For the purpose of later use, we write Qo, in Eq. (fT9i) . in the 
integral form in the 3D Euclidean momentum space, 

^ 2 lAI 2 



n (a,\A\) 



a 



+ 



AGs 4G D 



E f d 3 p E 
, J W 



ln[^±^U21n^ + ^ 2 



Pe 



Pe 



(21) 



The p E terms in the denominators in the second line of Eq. (l2~Ii) are inserted so as to drop 
an irrelevant infinite constant. 



B. Renormalized thermodynamic potential 

As mentioned in the previous subsection, Q , T = part of Q is ultraviolet divergent. To 



eliminate the divergences, we introduce the counter Lagrangian as derived in 



(see also 



fill): 



C c =~Z s a 2 -Z D \A\ 2 , (22) 

6 3 2 X 

Z s = —A- -a, Z D = —A- -a, (23) 

7l 2 4 TT 2 4 

where A is the 3D momentum cut-off and a is the arbitrary renormalization scale. For 
completeness, the derivation of Eq. (l23|) is given in Appendix B. 

Introducing the above counter Lagrangian, Q in Eq. (f2~Tj) turns out to be finite and the 
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renormalized Qq t becomes 



n 0r (a, |A|) 



4Gc 



- -a (T z + 



4G 



r )|A| 



(2tt) 3 



In 



P 2 eo + E 

Pi 



±2 



+ 2 In 



±2 



Pi 



2" " 

Pe 



Pi 



IAI 



(24) 



Note that the counterterms cancel the divergences and the integral becomes finite. Thus, 
after performing the renormalization, the renormalized thermodynamic potential Q r (= f2 0r + 
Qt) is finite and we carry out the numerical analyses on Q r . 

Before studying the qq and qq condensates, we rewrite Q r by using the following param- 
eters: 



2?r/ 1 



3 

-a 



A 



-7T 



4G 



1 

4° 



(25) 
(26) 



Using (Jq and A , we can write Q r at T = = /x as 



fi r (a,|A|) 



1 



T=0=At 
1 



CTo(T 2 _ _a |A| 



2tt 



TT 



+ —^ + —((7+^1) 



3tt 1 



A| 



(27) 



37T 37T 

Throughout in the following, we use the parameters <7o and Ao- By minimizing Eq.( l271) for 
A = 0, one can easily verify that, when cr > 0, cr corresponds to the qq condensate in 
vacuum. 



IV. QUARK- ANTIQUARK AND DIQUARK CONDENSATES 

We have obtained the thermodynamic potential in the previous section. Eqs. (1241) . fT25|) 
and (1261) tell us that this model has two free parameters (<To, A ). In this paper, we aim to 
study the system where the qq condensate always takes place in vacuum if A = 0, so we 
assume a to be positive (see, e.g. [jjj]). After fixing a , there remains a free parameter A 
and we introduce r through 

r = A /a . (28) 

There is no direct way of fixing the parameter r, and we analyze for different r's. 

Fig. [TJ plots the qq and qq condensates (normalized by do) at T = 0. In the case of 
r = —10 (panel (a)), we see that the qq condensate disappears at /i = 1.0cr , and there does 
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FIG. 1: a (circles) and A (triangles) as a function of chemical potential \i at T = 0. 

not arise the qq condensate. Through numerical analyses, we have found that this is the 
case for r < —6.3. On the other hand for r = — 1 (panel (b)), at /j, — 1.0a the qq condensate 
disappears and at the same time, the qq condensate arises. Similar results are obtained for 
r = and 1, where the transition densities are fi = 0.88o"o and \x = 0.4o"o, respectively. 
All the phase transitions in the panels (a)-(d) are of the first order. With increasing the 
ratio r, the qq condensate A becomes larger and eventually exceeds the qq condensate at 
/j, — 0. As seen from the panel (e) for r = 3/2, the qq condensate disappears and only the 
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r = -1, T = 0.1oo 




0.2 0.4 0.6 



1 1.2 1.4 1.6 




-1, T = 0.4ao 



0.6 
0.4 
0.2 

0. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



FIG. 2: The two gaps a (circles) and A (triangles) for (r, T) = (—1, O.lcro) and (—1, 0.4o"o) 



qq condensate exists for whole fi. More detailed analyses show that the qq condensate does 
not occur for r > 1.15. Thus the results are sensitive to the ratio r = A /cro- It should be 
emphasized that there is no region where the qq and qq condensates coexist. To clarify this 
fact, we show the close-up of the condensates near the phase transition point in Fig. [U(c). 

Now we turn to the T ^ case. We display the results for (r, T) = (—1, 0.1cr ) and 
(— 1, 0.4o"o) in Fig. El These panels show that the qq condensate for jj, — is 1.0cr at 
T = O.lcro and 0.92er at T = 0.4<Jo. The qq condensate for /i = l.Odo is 0.24o"o at T = O.lcro 
and at T = 0.4cr . Thus, as T increases, the qq and qq condensates decrease. Note 
that the qq condensate appears at T = O.lcro and does not appear at T = 0.4cro, which 
indicates that A disappears at high temperature. More detailed analysis shows that the qq 
condensate at // = 1.0cr disappears for T > 0.15<t (see Sec.V). It should be noted that, 
within our numerical accuracy, there is no region where the qq and qq condensates coexist, 
also at finite temperature. The results for other values of r are qualitatively the same. As 
the temperature increase, the condensates become smaller and completely disappear at the 
critical temperature. This is the signal of the phase transition from the condensate state to 
the normal state. 

From Fig. [2j we see that the phase transition for T = O.lcro is apparently of the first 
order and the case for T = 0.4cro, it is of the second order. Investigating the thermodynamic 
potential as a function of the qq condensate as in 12 ]. we have confirmed this fact. We will 
discuss the phase transition and its order in the next section in more detail. 
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FIG. 3: The phase diagram of the 3D GN model. 
V. THE PHASE DIAGRAM 

Through minimizing the thermodynamic potential, one obtains the phase diagram. In 
Fig. [3l we display the phase diagrams for various values of r. For r < —6.3, there appears 
the pure qq condensate phase at low temperature and density and no qq condensate phase 
appears. In the cases of r = —1, 0, the phase diagrams bear resemblance to that of QCD. 
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For r = — 1 with T = 0, the phase transition from the qq condensate phase to the qq 
condensate phase takes place at /i ~ 0.98o"o. For fi = 0, as T increases, the transition 
from the qq condensate phase to the normal phase takes place at T ~ 0.73<7o, which applies 
also for r = —10, 0, 1. On the other hand, in the cases of r = 3/2 and 10, the transition 
temperature for \i = are T ~ 1.08a and 10<t , respectively. 

As r increases, the region of qq condensate phase shrinks toward the \x axis and the 
region of the qq condensate phase increases toward the T axis. For r = 3/2 and 10, the qq 
condensate does not exist only and the qq condensate appears. Through numerical analysis, 
we have found that the qq condensate phase disappears completely at r ~ 1.15. 

The points (/x c , T c ) shown in the panel (a), (b) and (c) in Fig. [3] represent the critical 
points from the first order phase transition to the second order. The phase transition below 
the critical temperature T c is the first order and above T c is the second order. On the other 
hand in the panel (d), there is no critical point and the phase transition from qq condensate 
to the qq condensate is always the first order. The more detailed analysis tells us that the 
critical point disappears for r ~ 1, and it always appears between the qq phase and the 
normal phase. With respect to the qq condensate, the phase transition from the qq to the 
normal phase is always of the second order. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the qq and qq condensates in the 3D GN model with 2d spinor quarks, 
and obtained the phase diagram for various values of r = A /a . 

We have found that the behaviors of the qq and qq condensates at T = 0, in Fig. [TJ bear 
resemblance to that of the 4D NJL model With increasing r, the qq condensate becomes 
more dominant and the qq condensate disappears completely for r > 1.15. 

Fig. [1] and Fig. [2] show that, both for T = and T ^ 0, there is no region where the qq 
and qq condensates coexist. This is a characteristic feature in the 3D GN model which does 
not happen in the 4D NJL model. 

From Fig. [3] (b) and (c), we see that for r = —1, 0, there is a close resemblance be- 
tween the phase diagrams for r = — 1, and that of QCD. The qq condensate phase here 
corresponds to hadronic phase in QCD, and the qq condensate phase corresponds to color 
superconducting phase. In the case for r < —6.3, there does not appear the qq condensate 
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phase, and the phase diagram shows close similarity with the QCD phase diagram without 
color superconducting phase. However the diagrams for r = 1, 3/2 and 10 are very different 
from the QCD case. Especially when r = 3/2 and 10, there is no qq condensate and only 
the qq condensate exists. 

The circles in Fig. [3] indicate the critical points with respect to the qq phase transition 
from first order to the second order. Note that the phase transition from the qq condensate 
to the qq condensate is always of the first order. As seen from Fig. [1] and |2l when qq 
condensate arises, the qq condensate disappears rapidly. 

We have found that the phase structure drastically changes according to the value r = 
A /o"o. With increasing r, the qq condensate becomes larger. This means that the qq 
condensate becomes more dominant when A increases. On the other hand, A is related 
to the qq coupling constant Gd through Eq. (|26|) . and as A increases, Gd increases. In 
the same reason, when G$ increases (decreases), o"o becomes large (small), which causes r 
to decrease (increase). With these observations in mind, we can conclude that when Gd is 
small, the qq condensate is dominant over the qq condensate, and the qq condensate becomes 
larger when Gd increases. This is the same phenomenon as seen in the 4D NJL model. 

Finally it should be mentioned again that we have found the absence of the coexisting 
phase in the 3D GN model with 2d spinor quarks. As mentioned in Sec. I, the 3D GN model 
with 4d spinor quarks bears resemblance to the 4D NJL model. Then it is worth studying 
the phase diagram in the 3D GN model with 4d spinor representation to see whether the 
coexisting phase appears as in the 4D NJL model or does not appear as in the present case. 
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APPENDIX A: THE DERIVATION OF EQ.flHD 



The determinant of G 1 can be rewritten as 

<t + fij°)lfl e ier 2 Al s 



DetG" 1 = Det 



ier 2 A*l s {j-<r-fi>f)l f l c 



x Det | Mf-' + AiT )!/ ° | (A 1) 

y- a -^) 1} 

where l c is the unit matrix and e is the antisymmetric matrix in color (red and green) space, 
e rg = —€ gr = 1- For a 2 x 2 block matrix with matrices A, B, C and D, we have the identity 

I A B \ 

Det = Det (-CB + CAC- l D) . (A2) 

\C D J 

Replacing A, B, C and D with corresponding elements in the first line of Eq. (1All) . we get 

- a + /i7°)l/l c ieT 2 Al s 
ter 2 A*l s - a - iil )lfi 

= Det (-|A| 2 +p 2 -p 2 + a 2 - /i 2 -2a ^-^ fry + /ry fr)* 
= (pi -E 2 -ii 2 - |A| 2 - 2y/E 2 n 2 + a 2 \A\ 2 Y 

x (pi -E 2 -fi 2 - \A\ 2 + 2y/E 2 fj, 2 + a 2 |A| 2 ) 4 , (A3) 

which leads to the first two factors in Eq. (|14p . After calculating the second line of Eq. (IA1[) . 
we finally obtain Q(a, |A|) in Eq.p^. 



Det 



APPENDIX B: RENORMALIZATION IN VACUUM 



As mentioned in Sec. I, the standard O(N) 3D GN model is renormalizable in the leading 
1/N order. This is also the case for the thermodynamic potential Eq. (1161) in the present 
model with in the mean field approximation. 

Following the procedure as in jl3j] (see also 6|), we carry out the renormalization and 
obtain the renormalized thermodynamic potential. 

The divergent part of the potential is Qq in Eq. (I2ip and the divergent integral is written 

as 



n 



Odiv 



>,|A|) = - 



:<7 



! A 



|A| 2 A 



7T- 



7T 



(Bl) 
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It is to be noted that fiodiv is independent of /i. 

To eliminate this divergence, we introduce the counter Lagrangian Cc, Eq. (1221) . For 
determining Z$ in Cc-, we consider the one-loop radiative correction to the a propagator in 



the 2 + 1 Minkowski momentum space as calculated in 



13| 



D -W = w^rmfy (B2) 

H(p ) - -N f N c J (27r)3 k 2 {k _ p) 2 > (B3) 

where Nf and N c are the numbers of flavors and colors. II(p 2 ) shown above is ultraviolet 
divergent, and the counterterm Z$ is introduced so as to eliminate its divergent contribution. 
Here we employ the renormalization condition 

D a (p 2 ) = -2iG s , at p 2 -p 2 = -a 2 . (B4) 

This means that the counterterm Zs should satisfy the following relation 

Z s = -tU(p 2 )\ p ^ 2= _ a2 . (B5) 

Going into the Euclidean space, and restrict the integration region with a sphere of radius 
A. Then after evaluating integral Eq. flB3l) and using the condition Eq. flB5l) . we obtain 

Z S = -^A- -a. (B6) 

7T 4 

In the similar manner, one obtains the renormalization constant Ze>: 



Z D = - 2 K- -a. (B7) 

7T 4 



Introduction of the counter Lagrangian Cc eliminates f^odiv, and we obtain the renormalized 
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